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ABSTRACT 


Numerical Method is the important aspects in solving real world problems 
that are related to Mathematics, science, medicine, business etc. In this 
paper, We comparing the two methods by using the scilab 6.0.2 software 
coding to solve the iteration problem, which are Gauss Jacobi and Gauss 
Seidel methods of linear equations. 
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INTRODUCTION 

A Linear System of variables x 1 ,x 2 , . x n is linear equations 

of the form a 1 x 1 +a 2 x 2 + . CL n x n - b n . The value of 

a 1 ,a 2 . a n are any constant real/complex numbers. The 

constant is called the coefficient of x 1 and b is called the 
constant term of the equation.A system of linear equations 
( or linear system } is a finite collection of linear equation 
in same variables. 

a ll x l +a 12 x 2 + . a ln x n = bl 

a 21 x l + a 22 x 2 + . a 2n x n~ ^2 

Q"n\ x l+®'n2 x 2 + . ^nn x n~ bn 

A solution of a linear system is n- tuple (tq, v 2 , . v n ] of 

numbers that satisfied each linear equation when the 

values v 1 ,v 2 , . v n are substituted for x 1 ,x 2 , . x n 

respectively. 

There are two Iterative methods for the solving 
simultaneous equations. 

1. Gauss Jacobi Method 

2. Gauss Seidel Method 

It can be shown that the Gauss-Seidel method converges 
twice as fast as Jacobi method. 
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METHODS 

1. Gauss Jacobi Method 

Let us explain this method is the case of linear equations 

in three unknowns. 

Consider the system of linear equations, 
a 1 x+b 1 y+c 1 z= d ± 
a 2 x+b 2 y+c 2 z= d 2 
a 3 x:+b 3 y+c 3 z= d 3 

Forconverges |<Zil> \bi\ + |c x | 

\b 2 \> \a 2 \ + |c 2 | 

|c 3 |> |a 3 | + |b 3 | 

For iterative process 

x =^(di- Ky- c t z] 
y-~[d 2 - a 2 x- c 2 z ] 
z=z~[d 3 - a 3 x- b 3 z) 

c 3 

If x (0) ,y (0) ,z (( P are the initial values of x,y,z respectively, 

then 

y (1) =^-(d 2 - a 2 x m - c 2 z (0) ] 
zD)=—(d 3 - a 3 xb 3 z 

c 3 
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Again using these value x^.y^.z^ we get 
x^=^-{d 1 - b ± y^- c-lZ^) 
y m =^[d 2 -a 2 xW-c 2 z^ 
z (2) =—(d 3 - a 3 x (1) - b 3 z (1) ] 

c 3 

Proceeding in the same way , if the rth iterates are 
x (r) ,y (r) ,z (r) , the iteration scheme reduces to 

x (r+i)_J_h 1 yf r )- CjZ*^] 

y (r+1) =f(d 2 -a 2 x (r) - c 2 z (r) ] 

°2 

z ( r+ i)=—[d 3 - a 3 x^- b 3 z (r) ] 

c 3 

The procedure is continued till the convergence is assured 
(correct to required decimals]. 


2. Gauss-Seidel Method 

This is only a refinement of Gauss-Jacobi Method As 
before, 

x =^ d i- b i y- c i z 3 

y =j-{.d 2 -a 2 x- c 2 z] 

°2 

z=^-(d 3 -a 3 x-h 3 y] 

c 3 

We start with the initial values y(°),z(°) for y and z and get 
That is, 

b t y^- c^ 0 )] 

We use z (0 ^ for z and for x instead of x (0) as in the 

Jacobi’s Method, we get 

y (1) =i-(d 2 -a 2 y (1) -c 2 z (0) ) 

°2 


PROBLEM: Solve the following equation by Gauss 
Jacobi Method. 

10x-5y-2z=3; 4x-10y+3z=-3; x+6y+10z=-3. 
Solution: 

The below result is find out by the manual work: 


X y z 


1 

0.3 

0.3 

-0.3 

2 

0.39 

0.33 

-0.51 

3 

0.363 

0.303 

-0.537 

4 

0.3441 

0.2841 

-0.5181 

5 

0.33843 

0.2822 

-0.50487 

6 

0.340126 

0.283911 

-0.503163 

7 

0.3413229 

0.2851015 

-0.5043592 

8 

0.34167891 

0.2852214 

-0.50519319 

9 

0.341572062 

0.285113607 

-0.505300731 


By using Scilab coding the result shown below: 

■< finallgaussjacobi.sci (C:\Users\Admin\Documents\finallgaussjacobi.sci) - SciNotes 
File Edit Format Options Window Execute ? 


finallgaussjacobi.sci (x) | 

1 

cl c;clea r ;close: 


2 



3 

x0=0;y0=0;z0=0; 


4 

deff ('x=fl(y,z) ' 

, 'x=(3+5*y+2*z>/10’ ) 

5 

deff ('y=f2(x,z) » 

# ' y=(-3-4*x-3*z)/-10' ) 

6 

deff <’z=f3(x,y) ' 

, 'z=(-3-x-6*y)/10') 

7 

for i=l: 9 


8 

xl=fl{y0,zO) 


9 

yl=f2(x0,zO) 


10 

zl=f3(x0,y0) 


11 


%i) = %g\n\n\ty( %i) =-%g\n\n\tz(%i)=%g\n\n\n\n' , i, xl, i, yl f i, zl ) 

12 

x0=xl;y0=yl; 

zO-zl; 

13 

end 


14 

princf ('Thus we 

find chat - solution - converges■to -%g,%g- and %g' , x0, y0, zO ) 

15 




n 

X 

y 

z 

i 

0.3 

0.3 

-0.3 

2 

0.39 

0.33 

-0.51 

3 

0.363 

0.303 

-0.537 

4 

0.3441 

0.2841 

-0.5181 

5 

0.33843 

0.28221 

-0.50487 

6 

0.340131 

0.283911 

-0.503169 

7 

0.341322 

0.285102 

-0.50436 

8 

0.341679 

0.285221 

-0.505193 

9 

0.341572 

0.285114 

-0.5053 


The Actual Values are x=0.3415, y=0.2852, and 
z=—0.5053. 


Now, having known x (1) and y (1) use x ® for x and y (1) for 


y ,we get 

z ( 1 ) =z-[d 3 -a 3 x ci) - bsy 1 ) 

c 3 


In finding the values of the unknowns, we use the latest 
available values on the right hand side. If x^ r \y^ r \z^ are 
the rth iterates, then the iterates scheme will be 

x (r+l)_J_^ 1 _ fo 1 y( r ). CjZ^] 

y( r+ i)=—(d 2 - a 2 x( r+ V- c 2 z (r) ] 
b 2 

z ( r+ i)=—(d 3 - a 3 x (r+1 )- h 3 y (r+1) ] 

c 3 

This Process of iteration is continued until the 
convergence is assured. 

PROBLEM: Solve the following equation by Gauss- 
Seidel Method. 

10x-5y-2z=3; 4x-10y+3z=-3; x+6y+10z=-3. 
Solution: 


The result is find out by the manual work. 


N 

X 

y 

z 

1 

0.3 

0.42 

-0.582 

2 

0.3936 

0.28284 

-0.509064 

3 

0.3396072 

0.28312364 

-0.503834928 

4 

0.34079485 

0.28516746 

-0.50517996 

5 

0.3415547 

0.28506792 

-0.505196229 

6 

0.3414947 

0.285039 

-0.5051728 

7 

0.3414849 

0.28504212 

-0.5051737 


By using the Scilab coding the result is shown below: 

•*•7 finalgaussseidel.sci (C:\Users\Admin\Documents\finalgaussseidel.sci) - SciNotes 


File Edit Format Options Window Execute ? 


Li Ea 11 □ ffii , ® * h | It | i 

finalgaussseidel.sci (x) 


1 

2 

3 

4 

5 

6 
7 


-:; clear: close: 

xO=j;yO= ;z0=:; 

deff ( 'x=fl (y,z)' 
deff ('v=f2 (x.z)• 
dgff {'z=f3(x,y)* 
for i=::" 


8 

9 

10 
11 
12 

13 

14 


x0=fl(y0,z0) 
y0=f2{xO,zO) 
z0=f3 Jx0,y0) 

end 


O') 

- 10 ') 

’) 


,i,x0,i,y0 f i,z0i 

nverges t: and hg ,xO,yO,zO) 
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N 

X 

y 

z 

1 

0.3 

0.42 

-0.582 

2 

0.3936 

0.28284 

-0.509064 

3 

0.33960 

0.283124 

-0.503835 

4 

0.340795 

0.285167 

-0.50518 

5 

0.341548 

0.285065 

-0.505194 

6 

0.341494 

0.285039 

-0.505173 

7 

0.341485 

0.285042 

-0.505174 


The Actual value are x=0.3415, y=0.2852, and z= - 
0.5053. 

CONCLUSION 

There are different methods of solving of linear equation 
some are direct methods while some are iterative 
methods. In this paper , two Iteration methods of solving 
of linear equation have been presented where the Gauss- 
Seidel Method proved to be the best and effective in the 
sense that it converge very fast with Scilab Software. 
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